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Non-Fermi liquids from holography 

Hong Liu, John McGreevy and David Vegh 
Center for Theoretical Physics, Massachusetts Institute of Technology, Cambridge, MA 02139 

We report on a potentially new class of non-Fermi liquids in (2+l)-dimensions. They are identified 
via the response functions of composite fermionic operators in a class of strongly interacting quantum 
field theories at finite density, computed using the AdS/CFT correspondence. We find strong 
evidence of Fermi surfaces: gapless fermionic excitations at discrete shells in momentum space. 
The spectral weight exhibits novel phenomena, including particle-hole asymmetry, discrete scale 
invariance, and scaling behavior consistent with that of a critical Fermi surface postulated by Senthil. 



ON 
O 

o 



43 



> 

(N 

o 
o 



X 



I. INTRODUCTION 

The normal state of the high-Tp superconducting 
cuprates, and metals close to a quantum critical point are 
examples of non- Fermi liquids, which have sharp Fermi 
surfaces but whose low energy properties differ signifi- 
cantly from those predicted from Landau's Fermi liquid 
theory 0, 0, Eg. While Landau Fermi liquids are con- 
trolled by a free Fermi gas fixed point with (almost) no 
relevant perturbations [J, [H, [H, 0| j a proper theoretical 
framework characterizing non-Fermi liquid metals, which 
likely involves strong couplings, is lacking. In this paper 
we search for new universality classes of non-Fermi liq- 
uids using the AdS/CFT correspondence HUES]. 

According to the correspondence, any (quantum) grav- 
ity theory in a (d + l)-dimensional asymptotically anti- 
de Sitter (AdS^+i) spacetime is dual to a rf-dimensional 
quantum field theory "living at its boundary" . Through 
the AdS/CFT dictionary a gravity theory can (in prin- 
ciple) be used to obtain all physical observables of its 
boundary dual, like the physical spectrum and correla- 
tion functions. Compared to conventional approaches, 
the gravity approach offers some remarkable features 
which make it a valuable tool for discovering new strongly 
coupled phenomena: 

1 . At small curvature and low energies known gravity 
theories reduce to a universal sector: classical Ein- 
stein gravity plus matter fields. Through the du- 
ality, this limit typically translates into the strong- 
coupling and large- N limit of the boundary theory, 
where TV characterizes the number of species. Thus 
by working with Einstein gravity (plus various mat- 
ter fields) one can extract certain universal proper- 
ties of a large number of strongly coupled quantum 
field theories 1 . 

2. Highly dynamical, strong-coupling phenomena in 
the dual field theories can often be understood on 
the gravity side using simple geometric pictures. 
Familiar examples include confinement and chiral 
symmetry breaking in a non-Abelian gauge theory. 



3. Putting the boundary theory at finite temperature 
and finite density corresponds to putting a black 
hole in the bulk geometry. Questions about com- 
plicated many-body phenomena at strong coupling 
can be answered by solving linear wave equations 
in this black hole background. 

Consider a quantum field theory which contains 
fermions charged under a global U(l) symmetry 2 . When 
a finite U(l) charge density is introduced into such a the- 
ory, it is natural to ask whether the system possesses a 
Fermi surface and if yes, what are the low-energy excita- 
tions. One approach to these questions is to study spec- 
tral functions of a fermionic composite operator 3 . The 
presence of a Fermi surface may be revealed by the ap- 
pearance of gapless excitations of the operator at discrete 
shells in momentum space. 

This will be the approach taken here. We obtain 
the spectral functions for fermionic operators using the 
classical Einstein gravity. We will work in (2 + 1) di- 
mensions and leave other-dimensional theories for future 
study. While the string theory landscape in principle pro- 
vides a large number of AdS/CFT dual pairs, only a few 
explicit examples are understood in detail; in (2 + 1)- 
dimensions, these include the Af = 8 M2-theory and 
recently-discovered AB JM theories [ll|, HH ■ As em- 
phasized earlier, however, by working with the classical 
Einstein gravity we are extracting universal properties of 
a large number of boundary field theories, even if their 
explicit Hamiltonians arc not known. The spectral func- 
tions we find give strong indications of the presence of 
Fermi surfaces of some non-Fermi liquid. We find poles 
representing 'marginal' quasi-particles at discrete shells 
in momentum space, with scaling behavior different from 
that of a Landau quasi-particle. We also observe some 
other novel phenomena, including particlc-holc asymme- 
try and discrete scale invariance for a continuous range 
of momenta. 



1 The field theory origin of such universality is still rather myste- 
rious. 



The theory may also contain charged scalars, and both scalars 
and fermions may couple additionally to some non-Abelian gauge 
bosons. 

The spectral function of an operator is a measure of the density 
of states which couple to the operator. It is proportional to the 
imaginary part of the retarded function Gr of the operator in 
momentum space. 
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Our investigation was motivated in part by earlier work 
of Sung-Sik Lee [IJj], which initiated the study of spec- 
tral functions of fermionic operators using a gravity dual. 
Our results differ from those of we believe the dif- 
ference lies in the implementation of the real-time holo- 
graphic prescription [15|, [H, [TtJ ■ 

The plan of the paper is as follows. In the next sec- 
tion we set up the framework for calculating the spectral 
functions of a fermionic operator at finite density using 
the gravity description. In sec. IIHI we discuss properties 
of the spectral functions, including scaling behavior near 
a Fermi surface. We conclude in sec. IIVI with a discussion 
of the interpretation of the results and possible caveats. 



II. SET-UP OF THE CALCULATION 

Consider a three-dimensional relativistic conformal 
field theory (CFT) with a global U(l) symmetry that has 
a gravity dual. Such a system at finite charge density can 
be described by a charged black hole in four dimensional 
anti-de Sitter spacetime (AdS4) [lH, with the current J M 
in the CFT mapped to a U(l) gauge field Am in AdS. 
A fermionic operator O in the CFT with charge q and 
conformal dimension A is mapped to the gravity side to 
a spinor field ip with charge q and a mass 



mR 



(1) 



where R is the AdS curvature radius. The spectral func- 
tion of O at finite charge density can then be extracted 
by solving the Dirac equation for ip in the charged AdS 
black hole geometry. Which pairs of (q, A) arise depends 
on the specific dual CFT. However, since we are working 
with a universal sector common to many gravity theories, 
we will take the liberty of considering an arbitrary pair 
of (q, A), scanning many possible CFTs. 



with 

Q 2 M 



A = pL (l - — ) , n 



9fQ 
R 2 r 



(4) 



where r$ is the horizon radius determined by /(Vo) = 
0, and \i can be identified as the chemical potential of 
the boundary theory. For calculational purposes it is 
convenient to use dimcnsionlcss quantities. Consider the 
rescaling 



r->r r, [t,x)^> — (t,x), A ^—A , 
r R z 

M — > Mt-q, Q^Qr 2 a 
after which the metric becomes 

^ ee g MN dx M dx N = r 2 (-fdt 2 + dx 2 ) + ~, 
R H / 



with now the horizon at r = 1 and 



(5) 



(0) 



, Q 2 l + Q 2 . / 1" 
./=! + — ^ — , A = fi I 1 - - ] , f-i = gFQ 



The dimcnsionless temperature is given by 
T=±-(3-Q 2 ) . 



(7) 



(8) 



The zero-temperature limit is obtained by setting Q = 
•\/3. At zero temperature, near the horizon the metric (J6j> 
becomes AdS2 x M 2 with the curvature radius of AdS^ 
given by 



R 



(9) 



B. Dirac equation 



A. Black hole geometry 



To compute the spectral functions for O we need only 
the quadratic action of tp in the geometry d6j)— d9j) 



The action for a vector field Am coupled to AdS4 grav- 
ity can be written as 



S = 



1 

2^ 



d 4 x- 



r?2 
9f 



_G_ 

If 2 



(2) 



where g 2 F is an effective dimensionless gauge coupling 4 . 
The equations of motion following from (f 2l) a re solved by 
the geometry of a charged black hole [18, Il9i] 5 



ds 2 



R 2 



(-fdt 2 + dx 2 ) + 



R 2 dr 2 



(3) 



Sspinor = / d a+i x^g~i^T M V M 4> - mW>) (10) 



where 

4, = V M = d M + ^ afcA/ r ab - iqA M (11) 

and uiabM is the spin connection 6 . Note that the Dirac 
action (fTU|) depends on q only through 



fig = fiq = g F qQ 



(12) 



4 It is defined so that for a typical supergravity Lagrangian it is a 
constant of order O(l) 

5 For a generic embedding of l(2} into 4d Af = 2 supergravity, this 
solution can be lifted to an M-theory solution 12011 . 



6 We will use M and a, b to denote bulk spacetime and tangent 
space indices respectively, and fj,, v ■ ■ ■ to denote indices along 
the boundary directions, i.e. M = (r,fj,). All indices on Gamma 
matrices refer to tangent space ones. For notational convenience 
below we will take m to be defined in units of l/R, i.e. mR — > m. 
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i.e. through the combination of <7f9- This is expected; 
(i is the minimal amount energy needed to add a unit 
charge to the system, thus for a field of charge q, the 
effective chemical potential is given by [i q . Below, for 
notational simplicity, we will set gF = 1 and treat q as 
a free parameter, but one should keep in mind only the 
product of them is the relevant quantity. 

To analyze the Dirac equations following from (fTQ|) . it 
is convenient to use the following basis 



1 -1 ) ' 



7^ J ' W I V- 



(13) 

where are two-component spinors and 7^ are gamma 
matrices of the (2+l)-dimcnsional boundary theory. 
Writing 



V'i = {-gg rr )-ie- iuJt+lk * x ' 



(14) 



the Dirac equation becomes 



(15) 



by E3' 



(19) 



Equations ()15p can be further simplified by choosing 
the basis 7 = iai, 7 1 = ci, 7 2 = 03 and setting &2 = 0, 8 
after which one finds two sets of decoupled equations 

J— (d r T m^/g^P) y± — =F« (fti — u)z T , (20) 
J — (d r ± my/grT) z T = ±i(k 1 + u)y ± (21) 



where we have written <f>± 
ratios 

<-,+ — — , 



!J± 
z± 



Introducing the 



(22) 



and using (jl9|) . the retarded Green function can be 
written as 



Gr = lim e 

e->0 



-2m 



£+ 
e_ 



(23) 



where one should extract the finite terms in the limit. 

It is convenient to derive flow equations directly for £± 
as in [3. From (HOD, we find 



with 



K^r) = (-u(r),ki), 



9i 



-gtt 



-(w + /*,(!--)) 



(16) 

Note that since asr— >oo,it^cj + ^,j,w should corre- 
spond to the difference of the boundary theory frequency 
from fi q , i.e. u> = corresponds to the Fermi energy 

To extract the retarded Green function for O, we need 
to solve (fT5"j) with the in-falling boundary condition at 
the horizon and to identify the source and the ex- 
pectation value for O from the asymptotic behavior of tp 
near the boundary. Such an identification can be carried 
out from the prescription of [l6|, [TtJ , which amounts to 
identifying -0+ as the source and its canonical momen- 
tum in terms of r-slicing (which is essentially ip— ) as the 
expectation value. More explicitly, (j>± have the following 
asymptotic behavior near r — > 00, 



</>+ = Ar n 



Br~ 



(17) 



'— 3r£b = -2mV^7C± T (h Tu)± (h ± u)£ (24) 
The in-falling boundary condition at the horizon implies 

£fc|r=l = *. (25) 

With the boundary condition (|25[) . one can now inte- 
grate (|2~4"|) to r — > 00 to obtain the boundary correlation 
function directly. Below we will drop the subscript 1 on 

*!. 

At zero temperature, the a; — > limit of equa- 
tions ([20| (f2Tj) and (|24|) is singular, since gtt then has 
a double zero at the horizon. As we will see this has im- 
portant consequences for the behavior of Gr near u> = 0. 
Also, at w = the in-falling boundary conditions (j2"5|) do 
not apply and should be replaced by 



£±|r=l,w=0 



k 2 



6 



(26) 



Note that the 1 / y6 factor multiplying (U 9 in (f2"6"| has the 
same origin as the one appearing in ([9]). 



with 

= lfh^ A b = ^^D, k u = (-(u + u„), h) ■ 
2m- 1 ' 2m + 1 ^ V V ^ qh ' 

(18) 

Suppose the coefficients D (corresponding to expectation 
value) and A (corresponding to source) are related by 
D = SA, then the retarded Green function Gr is given 



7 Here we assume m > 0. For m < 0, one simply exchanges the 
roles of A and D. For m £ [0, both quantization procedures 
are allowed. Also note that the factor 7 in H19H comes from 
C*R ~ ({0,0^}) while in perturbing the boundary action we 
add -i J"(i 3 a;(^oC-l-(5^o) with (5 = O+7 . As noted in QJ, we 
must choose the sign of the overall gravity action to be consistent 
with unitarity. 

8 Since the system is rotationally symmetric, there is no loss of 
generality. 



III. PROPERTIES OF SPECTRAL FUNCTIONS 

A. General behavior 

We now describe the properties of Gr obtained by solv- 
ing (f2~4"|) . First note that by taking k — ► —k the equations 
for £± exchange with each other, leading to 

G 22 (uj,k)=G 11 (uj,-k) . (27) 

Similarly by taking q — > —q,LU — > — a; we hnd that 

Gu(«,fc;-?) = -^ a (-w,fc;?). (28) 

So it is enough to restrict to positive k and q. One can 
also check that as |w|, |fc| 3> both components reduce 
to those of the vacuum. When m = 0, by dividing the 
equation for £ + in ([24]) by wc find that = — 
which implies that 

G 22 (io,fc) = - 1 ? n = 0. (29) 

Gil (a/, fcj 

Combining (|2~7| and (f2"9"| wc thus conclude that at k = 0, 

Gn(w,fc = 0) = G 22 (cj,fc = 0) = i, m = 0. (30) 

Further study of G^ is possible by numerically solving 
. We will first consider T = and will mostly discuss 
the massless case. The mass dependence will be discussed 
briefly at the end. There are several consistency checks 
on our numerics. Firstly, ImGn and ImG22 are both 
positive, which is a requirement of unitarity since the di- 
agonal components are proportional to spectral densities. 
For a fixed large k 3> fi q , ImGn has a linearly-dispersing 
constant-height peak at lu + fi q « — k and ImG 22 has a 
peak at lu + [i q m k, while both components are roughly 
zero in the region lu + \i q <G (— k, k) (see figure Q] and [2]) . 
This recovers the behavior in the vacuum, which is given 
by PHI 

Gn=-J^±^r, G 22 = \\ k j^0\ (31) 
\j k + (uj + le) y k - (lu + it) 

with now the divergences at uj = ±fc smoothed out into 
finite size peaks. 

B. Fermi surface 

As one decreases k to [i q and smaller, the behavior of 
Gr deviates significantly from that of the vacuum. For 
dcfiniteness, let us now focus on q = 1 (with fi q = v3)- 
In this case the finite peak of Im G 22 in the large k region 
develops into a sharp quasi-particle-like peak near kp = 
0.918528499(1) (see figure^). The behavior of ImG 22 in 
the region of small kj_ = k~kp and lu can be summarized 
as follows: 




■- : w 
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FIG. 1: Spectral function ImG22(w) at k = 1.2 < fi q (left 
plot) and k = 3.0 > fi q (right plot) for m = and q = 
1 (fig = y/3). The function asymptotes to 1 as \w\ — > oo 
as in the vacuum (|31|l . Right plot: The onset of the finite 
peak at u> as 1.2 ~ k — fi q roughly corresponds to the location 
of divergence at uj = k in the vacuum (|3ip . The function is 
roughly zero between u G (— fc — fj, q , k — fi q ), as it is in vacuum. 
Left plot: The deviation from the vacuum behavior becomes 
significant. 




FIG. 2: 3d plots of ImGn(w,fc) and ImG 22 (^,fc) for m = 
and q — 1 (fJ, q = In the right plot the ridge at k ^> fj, q 

corresponds to the smoothed-out peaks at finite density of the 
divergence at u) = k in the vacuum. As one decreases k to a 
value kF ~ 0.92 < [i q , the ridge in ImG 2 2 develops into an 
(infinitely) sharp peak indicative of a Fermi surface. 



1. For k± < 0, we find a sharp quasi-particle-like peak 
in the region uj < and a small "bump" (with 
a broad maximum) in the region lu > (see fig- 
ure [3]). This appears to indicate that there is a 
quasi-particle-likc pole in the left quadrant of the 
lower-half complex w-plane. As k± — > 0_, both 
the peak and the maximum of the bump approach 
ui = 0, their heights go to infinity, and their widths 
go to zero. By carefully examining when the peak 
and the bump meet we are able to determine the 
accuracy of k F = 0.918528499(1) to 10th digit. 

2. For k± > 0, one does not see a sharp peak along 
real w-axis for either sign of lu. Instead one finds a 
"bump" (with a broad maximum) on the lu > side 
and a smaller bump on the lu < side. See figure [5] 
In the limit k± — > 0+ , both bumps approach lu = 
and their heights go into infinity. 

3. The quasi-particle-like peak and various bumps can 
also be studied by plotting Im G 22 (k,uj) as a func- 
tion of k for a given lu (see the left panel of figure [4] 
for a plot at uj = —0.001). In the limit lu — > 0_, the 
height of the peak goes to infinity with its width 
going to zero. At exactly lu = 0, however, the func- 
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tions Im Gn and Im G22 become identically zero for 
k > ^| = -j- (see the right panel of figure 2]). This 

behavior can be understood from (j26|) and (|24|) 
as follows. For k > ^| (at m = 0), the bound- 
ary conditions pp]) become real and since (|24[) are 
real equations, Im Gu(lu = 0, fc) are then identically 
zero in this region. Note that kp > 

Re[G 22 ], Im[G 22 ] 



Re[G 22 ], Im[G 22 ] 
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FIG. 3: ReG22(t<->) (blue) and ImG22(w) (orange) at k = 
0.90 < kp. In ImG22, at u < there is a quasi-particle- 
like peak; for ui > 0, there is a much smaller 'bump'. As k 
approaches /cf, the peak and the bump approach uj = and 
their heights approach infinity. The dashed lines are the real 
and imaginary parts of the fit function (|36p . Although the 
real part slightly deviates from the fit, there is a qualitative 
match. 



Im(G 



Im(Gii) 




FIG. 4: Left: Plots of ImGn(fc) (dashed line) and ImG 2 2(fc) 
as a function of k at u) = —0.001 (m = and q = 1). A 
sharp peak in ImG 2 2(fc) is clearly visible near kF ~ 0.9185. 
The height of the peak is finite. In the limit u — > 0_, the 
height of the peak goes to infinity and the location of the peak 
approaches kp from left. Right: Plots of ImGn(fc) (dashed 
line) and ImG22(&) as a function of k at u) — 0. For uj = 0, 
both functions become identically zero in the region k > ^| = 

Since kjp > l/v2 ; at lj = 0, ImG22 is identically zero 
around kF- 



Denoting the location of the maximum of the quasi- 
particle-like peak as w*(fcj_) we find that Lu*(kj_) scales 
with k± — > 0_ as 
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FIG. 5: ReG22(<j->) (blue) and ImG22(ci;) (orange) at k = 
0.925 > kF- One finds a "bump" at w > and a much smaller 
"bump" at u) < 0. As fcx approaches 0+, both bumps ap- 
proach lj = and their heights approach infinity. The dashed 
lines are the real and imaginary parts of the fit function ((39} . 
The fit is not so good for u < 0, though the qualitative trend 
matches. 



and the height of Im G22 at the maximum scales as 

ImG 22 (w,(fci),fci)~fcJ a , a =1.00 ±0.01 (33) 

(see figure |6|) . One finds exactly the same scaling be- 
havior also for the maxima of the other three "bumps" . 
This strongly suggests that in the limit of small k^_ and 
uj ImG22(^, has the following scaling form 

ImG 22 (A 2 w,Afci) = A- a ImG 2 2(w,jfe x ) (34) 

with the scaling exponents a, z given by 



2.09 ±0.01, 



a = 1.00 ±0.01 



(35) 



The scaling behavior (|32|) (|35|) suggests an underly- 
ing sharp Fermi surface with Fermi momentum kp. It 
is, however, not of the form corresponding to a Landau 
Fermi liquid which would have exponents z = a = 1. It 
is an example of the more gen eral scaling behavior dis- 
cussed recently by Senthil [23l . |24| for a critical Fermi 
surface occurring at a continuous metal-insulator tran- 
sition 9 . The system also has a rather curious particle- 
hole asymmetry; the quasi-particle-like peak at k± < 
morphs into a "bump" at k± > as the Fermi surface is 
crossed (compare feature 1 and 2 of previous page). The 



w*(/c_l) ~ k 



2.09 ±0.01 



(32) 



For early work, see [25, 26 



6 



fact that fi q > Uf suggests that the system has repulsive 
interactions. 

For a given uj ^ 0, Gu{k,ui) are non-singular for any 
value of k including kp, while for a given value of fc, 
Ga(k,uj) is continuous but non-smooth at u> = 0. 10 This 
non-smooth behavior at to = for momenta away from 
the Fermi momentum is puzzling and it would be nice to 
understand its physical interpretation better. From the 
gravity side, this is related to the aforementioned singular 
behavior of equations (f2"0|) - (j2l"j) and (|24|) near 10 — ► 0, as 
discussed around (|26|) . This can be further attributed to 
the existence of an AdS2 region © in the bulk geometry 
at zero temperature. 
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FIG. 6: Left plot: Dispersion relation aj*(fc) ~ k\ around 
kF- The dashed lines indicate the u- values of the maxima 
of the bumps. The solid line shows the dispersion of the 
quasi-particle peak. Right plot: Convergence of the z scaling 
exponent for k > kF (red) and k < kF (blue). The horizontal 
axis is the natural log of k±. 



Further support for the scaling behavior near the 
Fermi surface can be obtained by fitting the whole curves 
of G22(w) (rather than just the behavior near the max- 
ima) for different k± by a scaling function which is ana- 
lytic in the upper half u-plane 11 : 

1. For k± < 0, G22 can be fitted 12 for both signs of w, 
by (see figure [3]) 



G 22 (w, fcjj 



c (-kj_y 



log 



c 1 (~k±) z 



(36) 



where 7 0.34 and Co, ci are positive constants (in 
the scaling region). The above function has a pole 
in the lower half w-plane at 



-ci(-fcj_) 



(37) 



As k± — > 0_, lo c approaches to the origin of the 
complex plane along a straight line which has an 



From numerical calculation, it does appear that the functions 
become smoother for k 2> fJ-q- 

11 We would like to thank S. Sachdev for a discussion of possible 
subtleties. 

12 The function fits well not only along the real aj-axis, but also in 
the upper half plane. There are numerical instabilities at zero 
temperature in the lower half complex aj-planc and we have not 
been able to perform a direct fit there. 



angle 7 — n with respect to the positive real axis. 
Since Rew c gives the location u>*(k±) of the peak 
and —lmuj c gives the width T of the peak, for (1361) . 



tan 7 |o;*(£;_l)| . 



(38) 



This linear dependence of the width on w* 13 is rem- 
iniscent of the behavior in e.g. 0]. 

2. For k± > 0, we have not found a good global fit for 
both signs of u>. A reasonable fit for the imaginary 
part is 



G 2 2(w, kj_) 



a Q k A 



, M 



(39) 



where ao and ai are positive constants which take 
different values for u> < and ui > 14 . 



It is important to note that the functions ([36|) and (|39|) 
are only best numerical fits we could find and should not 
be taken too seriously as the "genuine" functions which 
describe the system. Both are consistent with the re- 
quirement from figure S] that ImG22(fc_i_) = at w = 0. 
The different fit functions for fcj^ < and k± > may 
reflect the "particle-hole asymmetry" discussed earlier. 
Also note that for a nonzero k± both (f3"6"|) and (f3"9"|) indi- 
cate a branch point singularity at to = 0, but have differ- 
ent k± — > limits. The behavior of ImG22(w) at exactly 
k = kp is not completely clear to us at the moment. 



Discrete scale invariance 



In the region k < where Im Gu(u 



-y^, ..iivn, iinujiv 1 " — 0, k) are 
nonzero (see figure [4]) , a new phenomenon occurs in the 
to ^ limit. One finds that lYnGu{uj,k) become os- 
cillatory with oscillatory peaks periodic in log \ ui\ with 
constant heights, see fig. [71 More explicitly we find 

G ii (LO,k) = G ii (uje n ^ k \k), neZ, u -» (40) 

where £(fc) is a (fc-dependent) positive constant. In other 
words, Ga is invariant under a discrete scaling in to. £(£;) 
appears to decrease with k and approaches a constant in 
the limit k — > 0. In the limit k — ► approaches in- 

finity. The height of the oscillatory peaks also decreases 
with k, approaching zero as k — > (where the whole func- 
tion approaches unity) and a finite constant as k — > 

It would be interesting to understand whether (|40[) is as- 
sociated with some kind of complex scaling exponents. 



13 Recall that for a Landau quasi-particle, r ~ a; 2 . 

14 Again we are handicapped by a numerical instability the lower 
half cj-plane which prevents a fit directly around the singularities 
in the lower half plane. 
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Below we will refer to the region k < as the oscil- 
latory region. Note that the oscillatory region appears 
to be the counterpart for fermions of the unstable region 
for a charged boson, where the corresponding bosonic 
modes have complex energies and want to condense. In 
the fermion case, the oscillatory region does not appear 
to indicate an instability, e.g. there is no singularity in 
the upper half of the complex w-planc. 



Re[G 22 ], Im[G 22 ] 
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FIG. 7: Both ReG 22 (^,fc = 0.5) (blue curve) and 
ImG22(i^,fc = 0.5) (orange curve) are periodic in logtj as 
lj — > 0. The period appears to be given by A(logw) ~ 

; a a . This formula was guessed based on the be- 

havior of the solution in the AdS2 region; the formula is con- 
firmed by the numerics. 



D. Finite temperature 

Turning on a small temperature T appears to smooth 
everything out. There is no longer a sharp Fermi surface, 
i.e. there no longer exists a sharp momentum at which 
ImG22 becomes singular for any real to and k. Going to 
the lower half w-plane, one finds that all the singularities 
arc a finite distance away from the real axis, with the 
closest distance given by T which happens at k w 0.90 
(sec figure [5]) 15 . This behavior is different from the Fermi 
liquid where the width is quadratic in temperature. Note 
that for a given small k± < 0, as one turns on the tem- 
perature, the corresponding quasi-particlc-likc pole in the 
complex w-plane appears to move down and to the right. 
It is also interesting to note that at finite T, there are now 
quasi-particle-like poles for momenta k > kp. Perhaps 
they are generated from the branch point at T = 0. 

At finite T, the functions Im G M become smooth at 
cu — and in the oscillatory region there are only a finite 
number of oscillations as the u> — > limit is approached. 
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FIG. 8: The complex omega plane for T = 4.13 x 10~ 4 : now 
the quasi-particle pole is finite distance below the real a>-axis. 
The dashed line indicates the trajectory of the pole between 
k = 0.87(left). . . 0.93(right). The closest distance to the real 
axis is equal to the temperature T (up to 1% accuracy). There 
is a numerical instability for Imui < —nT which can also be 
seen directly from the wave equation. We leave it for future 
work to explore this part of the lower half plane. Also shown 
is the density plot for ImG22(<^) at k = 0.90, where the cor- 
responding pole is closest to the real axis. 



E. Charge dependence 

When we increase (decrease) q to be greater (smaller) 
than 1, the Fermi momentum kp increases (decreases) 
with q approximately linearly. As q is further increased, 
new branches of Fermi surfaces appear. These features 
can be seen in figure [5J which gives the density plots of 
ImGn and ImG22 in the q — k plane at a fixed value of 
lu = -0.001. 

We have sampled the exponents z, a for a few other val- 
ues of q for the lowest branch of fermi surface in Im G22 , 
e.g. for q = 0.6, z ~ 5.32, a f» 1.00, and for q = 1.2, 
z w 1.53, a ?s 1.00. Compared to the values for q = 1 
described earlier, it then appears that z decreases rapidly 
with increasing q, while a = 1 is independent of q. Note 
that in [23| it was argued that z > a and z > 1. Thus it 
could be that z will asymptote to 1 for larger values of 
q. 16 We also find that the constant 7 in (|36| appears to 
decrease rapidly with q. Thus it seems likely that as q 
is increased, the non-Fermi liquid will become more like 
a Landau Fermi liquid. Given that kp increases with q, 
this is reminiscent of asymptotic freedom in high-density 
QCD. 

The q—k space in figurc[5]is separated into two regions 
by the (black) k = ^| line. In the region to right (sta- 
ble region), the locations of the quasi-particle lines (i.e. 
orange lines in figure [9]) stabilize in the limit uj —> and 
indicate locations of Fermi surfaces. The region to the 
left is the oscillatory region discussed earlier, where the 
log-periodic oscillatory behavior is reflected in a down- 
ward motion of the orange lines as |u>| is decreased; they 
seem to become infinitely dense in the limit u> — » 0. Also 



lj Similar results have also been obtained by Carlos Fuertes. We 
thank Carlos Fuertes and Subir Sachdev for communicating the 
results to us. 



16 At larger values of q the convergence of the exponents becomes 
slower; we leave this for future work. Also note that the value 
a = 1 is special according to the scaling theory of [23ll . 
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recall that in the oscillatory region, the heights of the 
peaks remain finite in the cu — > limit. 

As one decreases q, a Fermi surface line in figure [9] will 
intersect the line k = disappear into the oscillatory 
region, and lose its status as a Fermi surface. Thus the 
behavior of ImGu in the oscillatory region is strongly 
correlated with the sprouting of new branches of Fermi 
surface as q is varied (see figure [T0|) . 



Im(Gn) Im(G 22 ) 
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FIG. 9: Density plot of Im Gij (k, <?)o;=-o.ooi with k G [0,5] 
and q e [0, 7] at T = 2.76 x 10~ 6 . A negligible temperature 
was turned on in order to increase the speed of the compu- 
tation. The results were not affected by this. The orange 
lines are locations where the functions become very large. 
Also note that the width of the peaks in ImG(fc) decreases 
quickly as one moves towards larger charge. The black line is 
k — Hq/^/E to the left of which is the oscillatory region. 



Im(Gn), Im(G22) Im(G rl ) 




FIG. 10: (i) ImGy(/e)a, = o with q = 1.6. Near k = /i q /V&, 
a bump is seen in lmG22{k). At slightly higher value of 
the charge, the Fermi surface crosses the boundary of the 
oscillatory region and this bump becomes a peak. (ii) 
ImGii(<7, k) u =o along the fi q = V&k line. The spacing be- 
tween the peaks is constant, Ak ~ \/7/2(±l%). 

We have also studied other values of the mass (with 
m < i) at q — 1, and find the Fermi momentum kp de- 
creases linearly with increasing mass. At finite mass, the 

2 

oscillatory region is now given by k 2 +m 2 < (see ([26])). 
We find that for q — 1, the Fermi surface disappears into 
the oscillatory region at roughly m ps 0.4. For m > 
the oscillatory region disappears; we expect the Fermi 



surface will also disappear beyond this value if not be- 
fore. 



F. Summary 

To conclude this section, let us summarize the main 
features of the spectral functions which we have observed. 
The system we study is conformally invariant at zero den- 
sity. Turning on a finite charge density breaks Lorentz 
and scale invariance. The energy scale of the problem is 
controlled by the chemical potential [i which for a charge 
q particle becomes n q = [iq. At q = 1 the spectral func- 
tions also exhibit two other interesting scales. The first is 
the Fermi momentum kp < fi q around which we observe 
a quasi-particlc-like peak which suggests an underlying 
Fermi surface. The scaling behavior and the particle- 
hole asymmetry around the Fermi surface indicate that 
this is a non- Fermi liquid. The other scale is ks = 
which lies below kp. We find that for k < ks, the spec- 
tral functions have log-periodic oscillatory behavior near 
lu = 0, which indicates some underlying discrete scale in- 
variance. At larger values of q new scales corresponding 
to more than branches of Fermi surface also appear. 

It is important to emphasize that the scaling behavior 
observed here is not related to the scale invariance of the 
vacuum, which is broken by the nonzero charge density. 
It is emergent, arising as a consequence of the collective 
behavior of many particles. Note that both the scaling 
behavior around the Fermi surface and the discrete scale 
invariance involve the small uj limit, which on the gravity 
side can be attributed to the AdS2 region in the near 
horizon geometry of the black hole at zero temperature. 
It may be possible that this emergent scaling behavior 
can be understood from the AdS2 region 17 . 



IV. DISCUSSION 

Finally, we discuss some caveats and possible interpre- 
tation of the results. While the black hole geometry © 
is by itself thermodynamically and perturbatively stable, 
when it is embedded into a specific gravity theory, new 
instabilities may occur. For example, it is possible for 
charged (bulk) scalars to condense, spontaneously break- 
ing the U(l) symmetry [2|| |2{|. This happens if the bulk 
spectrum includes a charged scalar of sufficiently large 
charge or sufficiently small mass [3(| ■ The boundary the- 
ories considered in [301 ] all contain such scalars including 
the Af = 8 M2 brane theory and ABJM theory. It would 
be very interesting to understand how the condensate af- 
fects the Fermi surfaces and scaling behavior observed 
here, and how generic the existence of such scalars is. 



Work in progress with T. Faulkner [27|] , 
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The black hole solution ([6]) has a finite entropy at zero 
temperature, and thus describes an ensemble of an expo- 
nentially large number of states. Given that the solution 
is not supersymmetric, beyond the gravity approxima- 
tion likely these states are energetically closely-spaced, 
rather than exactly degenerate. This "frustrated prop- 
erty" is shared by many known models of spin liquids. 
The behavior described above then reflects average be- 
havior of a large number of states rather than that of a 
single ground state. 

A rather mysterious feature of our results is that differ- 
ent probe operators appear to find different Fermi surface 
structure which depends on (and only on) their charges 
and operator dimensions. One possible explanation is 
as follows 18 . Let us look at the OPE of e.g. the first 
component 0\ of a fermionic operator 0, which has the 
schematic form 

0^0,(0)--^ + ^-^ + ..., (41) 

where J° is the zero component of global U(l) current. 
(|4Tj) implies that the density for 0\ can be written as 

n 0l = (O^OM) ~ • • * - 

(42) 

where e should be considered as a short-distance cutoff. 
The first term on RHS of (f4"2"| is the standard piece due 
to vacuum fluctuations, which can be subtracted. In a 
state of finite density, the second term induces a density 
for 0i, which is proportional to q, the background charge 
density and depends on the UV cutoff through the dimen- 
sion of the operator. Since the Fermi surface involves 
modes with wavelengths of order kp , which is paramet- 
rically distinct from the short wavelength of the modes 
contributing the UV divergence, we expect the fermionic 
density which is responsible for the Fermi surface should 
be insensitive to the UV cutoff. The background charge 
density ( J°) oc N 2 induces a nonzero charge density of 
order 0(N°) for each charged operator, which depends 
on its charge and conformal dimension. The induced den- 
sity increases with its charge, which is consistent with our 
empirical observation that kp increases with the charge 
of the bulk particle. 

In the large N limit, the effective interactions of 
with itself and other gauge invariant operators are all 
suppressed by 1/N. This may lead one to conclude that 
the effective theory for should be a free Fermi gas, 
which would contradict our observed scaling behavior 
near kp. Note, however, the effective dynamics of 
can be different from that of a free fermion, again due 
to large N effects. To see this, let us look at the current 
density from 0, which can be schematically written as 



= 0y0. The fluctuations of which can be read 
from its connected two point functions, are suppressed by 
\/N 2 . Thus in the large N limit, does not fluctuate. 
One can try to model this by coupling a free fermion to 
a gauge field, which acts as a Lagrange multiplier sup- 
pressing the fluctuations of the associated current. When 
coupling a Fermi liquid to a dynamic al gaug e field, it is 
well known (sec e.g. 

that long-range magnetic interactions result in a non- 
Fermi liquid, which appears to be consistent with our 
picture. It would be desirable to make this argument 
more precise. Note that the particle-hole asymmetry is 
not seen in previously known models. The above sugges- 
tion does not preclude the existence of some fundamen- 
tal non-Fermi liquid structure from which the behavior 
of probe fermionic operators could be derived. 

The fact that the induced charge density for each probe 
operator is of order O(l) also implies that their contri- 
butions to the transport of the system are not visible 
at leading order in the large N expansion. Indeed, to 
leading order in N none of the observables like specific 
heat, conductivity, entanglement entropy can depend on 
the charge or dimension of probe spinor fields. However, 
if there exists a fundamental non-Fermi liquid structure, 
some effects might still be visible at leading order. We 
will leave this for future study. 

Finally, as indicated earlier, the near horizon AdS2 re- 
gion appears to play a role for the appearance of both 
the log-periodic behavior in the oscillatory region and 
the Fermi-surfaces. Clearly it would be interesting to 
have a better understanding of the CFT interpretation 
of the AdS2 region. 
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